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Abstract 

(N. 

■ We work out the semileptonic A^ — t- A£^£~ transition in standard as well as 

■ different supersymmetric models. In particular, considering the parametrization of 
CN ■ 

CS| \ the matrix elements entered the low energy effective Hamiltonian in terms of form 

OO \ factors in full QCD, we calculate the amplitude and differential decay rate responsible 
O. 

CSl ■ for this decay channel in supersymmetric models. We then use the form factors 
calculated via light cone QCD sum rules in full theory to analyze the differential 

, , branching ratio and lepton forward- backward asymmetry of this decay channel in 

^ I different supersymmetric models and compare the obtained results with those of the 



standard model. We also discuss how the results of different supersymmetric models 
deviate from the standard model predictions and which SUSY scenarios are favored. 

PACS number(s): 12.60.-i, 12.60.Jv, 13.30.-a, 13.30. Ce, 14.20.Mr 



*e-mail: kazizi@dogus.edu.tr 

^e-mail: sehban@istanbul.edu.tr 

■l-e-mail: a.t.olgun@gmail.coni 

^e-mail: z.tavukoglu@gmail.coni 



1 Introduction 



Recently, there has been an important progress on the course of searching for Higgs boson 
as a missing ingredient of the standard model (SM). The ATLAS and CMS Collaborations 
at CERN reported their observation on a Higgs-like particle with a statistical significance 
of 5(7 [1]. Now, it is searched whether this Higgs-like boson is the standard or non-standard 
Higgs particle. The supersymmetry (SUSY) has been the most popular paradigm for new 
physics (NP) scenarios in the last decades. The recent progresses have stimulated the 
theoretical works dedicated to the study of how a relatively heavy Higgs constrains the 
parameters of SUSY (for a discussion see for instance [2]). On the other hand, with these 
developments, we hope that we will have an experimental progress in searching for SUSY 
particles both directly by increasing the center of mass energy and indirectly by studying 
the flavor changing neutral current (FCNC) transitions. 

In the present work, we theoretically analyze the semileptonic FCNC decay of the At, — > 
A£'^£~ in existing related different supersymmetric models. In principle, the SUSY particles 
can contribute to such loop level transitions. Hence, we look for the effect of superparticles 
in this channel via calculating some related observables like differential branching ratio 
and lepton forward-backward asymmetry (FBA). Due to the specific features, there are 
different SUSY scenarios such as SUSY I, SUSY II, SUSY HI and SUSY SO(IO) [3-6]. 
In these models, the Wilson coefficients receive contributions from neutral Higgs bosons 
(NHBs) that are proportional to tan^/3, where tan/3 has been defined as the ratio of the 
vacuum expectation values of two neutral Higgs bosons {h^, A^). According to the tanP and 
an extra parameter fi with dimension of mass corresponding to mass term mixing of two 
Higgs doublets, the different SUSY models are categorized. In SUSY I, the /i takes negative 
vahie, some of the Wilson coefficients change their signs and the contributions of NHBs have 
been neglected. In SUSY II, the tan/S takes large value while masses of the superparticles 
are small in order of a few hundred GeV. In SUSY HI, the tan/3 is large and the masses 
of the superparticles are relatively large up to 450 GeV or more. In SUSY SO (10) model, 
the imaginary parts of the Wilson coefficients are large and the NHBs contributions are 
considered. 

In the last year, the CDF Collaboration at Fermilab [7] has reported the first observation 
on the baryonic FCNC transition of A° — >■ A//+/i~ with 24 signal events and statistical of 5.8 
deviations. They have measured a branching ratio of [1.73 ±0A2{stat) ±0.55{syst)] x 10~^. 
This decay channel is in the focus of different experiments fike LHCb at CERN [8] . Hence, 
theoretical and phenomenological predictions on the observables defining this channel can 
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help us in the course of searching indirectly for SUSY particles in this stage. Comparison 
of different theoretical results with experimental data may help us get useful informations 
about the existence of the SUSY particles. Note that the rare — >■ kt't' transition was 
analyzed in the same frameworks in [5] using only two form factors calculated via heavy 
quark effective theory (HQET). In this work, we generahze those calculations to include all 
form factors in full theory. 

In the next section, introducing the effective Hamiltonian both in the SM and SUSY 
models, we calculate the amplitude of the decay under consideration in terms of twelve form 
factors enrolled to the transition matrix elements. In section 3, we calculate the formula for 
the differential decay rate in SUSY and numerically analyze it together with the branching 
ratio and lepton FBA. We also compare the obtained results on the considered physical 
quantities in different SUSY models with those obtained from the SM. The last section 
encompasses our concluding remarks. 



2 The effective Hamiltonian and transition matrix el- 
ements 

In the SM, the A;, — >■ kt^i" transition goes with the h — >■ st't' at quark level whose 
effective Hamiltonian is given by 



eff 



GFOlemVtbVt 



ts 



2^277 



(2.1) 



where Gp is the Fermi coupling constant, «em is the fine structure constant at Z mass scale, 
Vij arc elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix; and the Cj-^'^ , C^^^^ and 
Cio are the Wilson coefficients. Considering the contributions of the new operators coming 
from the new interactions induced by the NHBs exchanged diagrams, the supersymmetric 
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effective Hamiltonian can be written as 



^SUSY 



ts 



- 2mbC'/f\sta^,q''{l + 75)^^7^^ - 2mhC!f" ^sia^^q^il - 75)617' 

q^ 

+ Cq,s(1 + 75)6£75^ + C^,s(1- 75)6175^ 



(2.2) 



where the new Wilson coefficients, Cq^ and Cq^ exist in the all considered SUSY models, 
while the primed coefficients only appear in SUSY SO(IO) scenario. 

The amplitude is obtained by sandwiching the new effective Hamiltonian between the 
initial and final baryonic states, i.e.. 



^5^5^'" = (a(pa) I nl^sY I Mpa.)), 



eff 



(2.3) 



where pa^ and p\ are momenta of the Af, and A baryons, respectively. To proceed, we need 
to calculate the following matrix elements parametrized in terms of twelve form factors in 
full theory: 



(A(pa) I S7/.(1 - 75)6 I Ab{PA,)) = ma(pa) 



7M/i(9')+^^x^.gV2(g') + rt(g') 



- 7M75^i(g^) - i(^nylb(tg2{q^) - q^l^g^iq^ 



(2.4) 



(A(pa) I S7^(l + 75)6 I Afe(pAj) = 'Ua(pa) 



7/./l(?')+^^M^?72(g') + rt(g') 



(2.5) 



(A(pa) I sia^^q^'^i + 75)6 | J^hiph^)) = '"a(pa) 



l,fl{q') + ^^^.gV^■(?') + q'fliq 



^ f_^r^2N , u.r/ 2^ 



+ 7^75^^ (g') + i(7,.ulbq"g'2{q') + q^lbQ^iq') 



(2.6) 



3 



(A(Pa) I siai^^q^il - -f5)b \ A6(paJ) = «a(pa) 



1/ 2\ ■ u T/ 2\ LL Ti 2\ 

7^75^1 {q ) - i(y^iulbq 92{Q )- Q 7593 {q ) 



(2.7) 



(A(pa) I s(l + 75)6 I A6(paJ) = —uk{pk 



rub 



- ilbQiiq^) - n'^(ynvlbq''g2{q^) - q^lbQziq^) 



(2.8) 



and, 



(A(pa) I s(l - 75)6 I A6(paJ) = —Mpa) 



nib 



4h{q^)+iq'<J^.uq''f2{q^) + q^h{q^) 



+ ilb9i{q^) + iq^(^nulbq''92{q^) + q^l59z{q^) 



(2.9) 



where is the transformed momentum squared; and the Ma^ and U\ are spinors of the initial 

(T) (T) 

and final baryons. In the meantime, the and with i — 1,2 and 3 are transition 
form factors. 

Using the above transition matrix elements in terms of form factors, we find the super- 
symmetric amplitude as 



M 



Ai^Ae+e- _ GFOiemytb^t 



ts 



SUSY 



2A/27r 



+ 
+ 

+ 



UA{PA)ilMiR + l3iL] + ta^.q''[A2R + B2L] + q^lAsR + -B3^])«a,(paJJ {h^i) 
UAiPA){li.[DiR + £iL] + ia^,q''[D2R + EiL] + q^iV^R + E3L\)uaSpaS\ Wld) 
MpaMQiR + niL] + tq^a^,q''[g2R + H^L] + q^iG^R + n3L])uA,{pA,)] (ii) 

'MPA)U[iCiR + S,L]+iq''a^,q''[IC2R + S2L]+q''[K:sR + SsL])uA,{pA,)] (^75^)| , 

(2.10) 



where R = (1 + 75)72 and L = (1 — 75)72 and the calligraphic coefficients are found as 



A, = hC|^f^-g^C|ff--2m,- 



fIC'7"^ + glC'r"" , M = A(l ^ 2), A, = A (1 ^ 3) , 



flC'r"^ - , i32 = i3i (1 ^ 2) , i33 = B, (1^3), 



T^i — /iC'io ~ fl'iC'io, 



V2 = Pi (1 ^ 2) , P3 = 2^1 (1^3), 



^1 — /iC'io + 5'iC'io, 



£2 ^£i{1^2), Ss^ Si (1^3), 



nib 



Hi = — 
rrib 



1 r 



5i 



rub 



/iC+ +^iCq^], 52 = 5i (1 ^ 2) , 53 = 5i (1^3), 



(2.11) 



with 



eff- 



Clo — C'lo + C'lO) 



10 



Cio ~ C'lo ' 



(2.12) 



3 Differential decay rate, branching fraction and FBA 



3.1 The differential decay rate 



In this part, we calculate the differential decay rate for the decay channel under consider- 
ation. Using the aforementioned amplitude, we find the supersymmetric differential decay 
rate in terms of form factors in full theory as: 



%{S)+Tl{s)z + T2{s)z' 



(3.13) 



where A = A(l, r,s) = (1 — r — s)^ — 4rs is the usual triangle function with s — (}^/m\^ 



r — mj^/mi and v — Jl 



Ami 



Here also z — cos 6 with 6 is the angle between momenta 



of the lepton Z+ and the in the center of mass of leptons. The calligraphic, %{s), 7i(s) 
and 72(5) functions are obtained as: 



%{s) = Z2m]mi^s{l + r - s) ( + l^sl' ) 
+ Um]m\^{l - r - s) Re V^E^ + Pg^:*] 



-DIE, 



+ (1 - r + s)Re 



+ 2>2mli2ml + mis) {{l-r + s)mK^^Re AIA2 + BIB2 



- mA,(l - r - s) Re AIB2 + A^B^ - 2^{rq A\Bi 



+ m\,s Re 



AIB2 



]) 
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+ 8mi^ |4mKl + r - s) + m\ [(1 - rf - s'] | ( I A|' + \Br\^ ) 



A + (1 + r - s)s] + m\s [(1 - r)' 1-421' + 



- 8mi < 4m^(l + r - s) - ml 



+ Smi^s^;^! - 8mA,sVr Re P;£:2] + 4(1 - r + Re P*P2 + <?r^2 

- 4(1 - r - s) Re [p*f2 + I^2*^i] + ^a, [(1 - r)^ - s^j [p^]' + |^2| 



8m^^ 4m^ 



(l-r)^-s(H-r) ReVllCi-^SlSi 
+ (4m,^ - mi^5) [(1 - rf s{l + r)] ( |^,|^ + ) 



+ Am\^\^s^{Am^ — m\^s) Re 



- 8mi^s|2v^(4m^2 - mi^s) (1 - r + s) Re [e?*^3 + ^t^s 

+ 4m^^/f(l - r + s)Re r>*/C3 + £:i*53 + P^/Ci + f:*^!] 
+ 4m^(l - r - s)Re VIS3 + £i*/C3 + P^^i + fg/Ci] 
+ 2(1 - r - s)(4m| - ml s) Re k*?^3 + ^1^3 



(l_^)2_5(l + ^)j(^|/C,|2 + |5i|2)| 



- 32mi^v^s|2m^Re + S^JCi + (4m| - mlj) Re ^^^T^iJ 

+ 8m^^s2|4^j^g + 2mA,V^(l - r + s)Re[/Ct/C3 + 5*^3] 

+ 2mA, (1 - r - s)Re [/C*53 + 8*^3 

- {^m\ - m\s){\ + r - J) ( |^3|' + l^al' ) 



Amell + r — s)Re 



P*;C3 + ^3*^53 



8m^-\/rRe 



©353 + S^IC-s 



+ Smlj-^ {l + r-s)( I/C3I' + l^al' ) + 4x/fRe 



/<^3'^3 



Ti{s) = -^2mi^rrnV\v{l-r)Re[A\Gr + Bl'H^ 

+ 2mA, {BIV2 - B;Vi + A^Si - AIS2^ 

+ 2mA^meRe (^AlHs + BlQs - A^Ui - B^Gi^ 

+ 2>2m\^s v\fx\^K^{l - r)Re(A*2V2 - B^E-^ 

+ x/f i?e [a^Vi + AIV2 - BlSi - BIE^ 

- ^nuRe (aIGs + BlKs + A\Qx + B^U^ 

+ 'i2m\m(,^vs^Re [a^Q^ + BlU^ , 



(3.15) 



r2(s) = -8mi^t;2A(|^i|' + |i3i|' + |I?i|' + |^i|') 



(3.16) 



In order to obtain the differential decay rate only in terms of s, we fulfill integrate 
Eq.(3.13) over z in the interval [—1, 1]. As a result, we get 



3.2 The differential branching ratio 



%{s) + \r2{s) 



(3.17) 



Using the differential decay rate, in this subsection, we numerically analyze the differential 
branching ratio and calculate the values of the branching ratios at different lepton channels. 
For this aim, we need sum inputs which we would like to present them here. In Table 1, we 
present the masses [9] as well as the hfetime of the initial baryon [9] , some constants and 
elements of the CKM matrix. 

The main inputs in our calculations are form factors. These form factors are calculated 
via light cone QCD sum rules in full theory in [10]. The fit function for the form factors 
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O T J. T» J 

borne Input Farameters 


Values 




0.00051 GeV 




0.1056 GeV 


ITLr 


1.776 GeV 




4.8 GeV 




5.620 Ge\/ 


m\ 


1.1156 GeF 




1.425 X 10-12 s 


h 


6.582 X IQ-^^GeVs 


Gf 


1.1 X 10 GeK 




1/137 




0.041 



Table 1: The values of some input parameters used in the analysis, 
/i, /2, /a, 9i, 92, 93, fL fl, 92 and gj is given as [10]: 

frm9rii')] = . \, + . ^ , (3.18) 

q 



1 - 




where the fit parameters a, b and m^.^ as well as the values of the related form factors at 
g2 = in full theory are given in Table 2. Furthermore, the fit function of the form factors 
fl and gf is given by [10]: 

flil'Mil')] = . \ X - . " , (3.19) 

1- ' 




where, the parameters c, rnj^^ and m^f^ as well as the values of the corresponding form 
factors at = arc presented in Table 3. 

In our numerical analysis, it is important to emphasize that the Wilson coefficient Cg^^^ 
has been taken to contain also the long distance (LD) effects coming from the charmonium 
resonances. These effects are parameterized using the Breit-Weigner ansatz as [11-13]: 

no ^ %C'-t., r>^r:^-- (3.20, 

where, G*-°^ = 0.362 and Ki are the phenomenological factors. Here, my^ and Fv; are 
the masses and decay rates of the vector charmonia, respectively. In the present work. 



9 



a b m^jj form factors at = 



fl 


-0.046 


0.368 


39.10 


0.322 ±0.112 






0.0046 


-0.017 


26.37 


-0.011 ± 0.004 




h 


0.006 


-0.021 


22.99 


-0.015 ±0.005 




9i 


-0.220 


0.538 


48.70 


0.318 ±0.110 




92 


0.005 


-0.018 


26.93 


-0.013 ±0.004 




■ ga 


0.035 


-0.050 


24.26 


-0.014 ±0.005 




fl 


-0.131 


0.426 


45.70 


0.295 ±0.105 




fl 


-0.046 


0.102 


28.31 


0.056 ±0.018 




9l 


-0.369 


0.664 


59.37 


0.294 ±0.105 




T 

■niB 













Table 2: The parameters in the fit function of the form factors /i, /2, /s, gi, g2, 93, /l^, /J, 
g2 and gj as well as their values a,t — in full theory [10]. 

c mj^^. m"j:jf form factors at = 

/f -1.191 23.81 59.96 0±0.0 
gl -0.653 24.15 48.52 ± 0.0 

Table 3: The parameters in the fit function of the form factors fi and gf as well as their 
values at = in full theory [10]. 

we only take into account the two lowest resonances that are J/'^{ls) and '4>{'2.s). The 
phenomenological factors have also been chosen as = 1 and ^2 — 2. The masses, 
branching fractions and total decay widths related to the considered resonances are given 
in Table 4. 











J/V^(is) 


3.096 


5.55 X 10-6 


92.9 X 10-6 







Tabic 4: The values of masses, branching fractions and total decay widths related to the 
resonances J/ip{ls) and ip{2s) [9]. 

Considering the above mentioned resonances from J/ip family, we divide the allowed 
physical regions into the following three regions in the case of the electron and muon as 
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final leptons: 

Region I ; Amf < < {mj/^^is) — 0.02)^, 

Region II ; {mj,^^u) + 0.02)2 < ^2 < {jn^^^s) - 0.02)2, 

Region III ; {m^{2s) + 0.02)^ < < {niAt, — "f^Af- 

In the case of r, we have the following two regions: 

Region I ; Am^ < q^ < {m^{2s) — 0.02)^, 

Region II ; (m^(2s) + 0.02)^ < q^ < {niA^, - m^)^. 



Finally, we would like to present the numerical values of the Wilson coefficients used in 

mimcrical calculations in Tabic 5. 



Coefficient 


SM 


SUSY I 


SUSY II 


SUSY III 


SUSY SO(10)(^ = -1000) 




-0.313 


+0.376 


+0.376 


-0.376 


-0.219 




4.334 


4.767 


4.767 


4.767 


4.275 




-4.669 


-3.735 


-3.735 


-3.735 


-4.732 










6.5(16.5) 


1.2(4.5) 


0.106 + Oi(l. 775 + 0.002i) 










-6.5(-16.5) 


-1.2(-4.5) 


-0.107 + Oi(-l. 797 - 0.002i) 













0.039 + 0.038i 
















0.011 + 0.072^ 








m ^ ^ 




-0.075 - 0.67i 
















-0.247 + 0.242i(-4.148 + 4.074i) 















Table 5: The Wilson coefficients used in numerical calculations [3,4,6, 14]. In the values 
containing parentheses, the values inside the parentheses stand for the r lepton, while the 

values outside belong to the e and /x cases. The other values (without parentheses) refer to 
all leptons. 



Having given all the inputs, we now present the dependence of the differential branching 
ratio on s for the e, fi and r leptons in the SM and different SUSY scenarios in Figures 
1 and 2. From these figures which are plotted considering the central values of the form 
factors, we see that 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 



s s 

Figure 1: The dependence of the differential branching ratio on s for the Af, — )■ Ae+e^ and 
Afe — )■ Afi'^fi~ transitions in SM and different SUSY models using the central values of form 
factors. 



30 




0.40 0.45 0.50 0.55 0.60 0.40 0.45 0.50 0.55 0.60 



s i 

Figure 2: The dependence of the differential branching ratio on s for the Aj, — )■ Ar+r^ 
transition in SM and different SUSY models using the central values of form factors. 

• in all lepton channels, the predictions of the SUSY II deviate maximally from those 
of the SM and other considered SUSY models. In the case of r, this deviation reaches 
to approximately one order of magnitude. 

• As far as the e and fi are concerned, the results obtained via the SUSY I have also 
considerable deviation from the predictions of the other models. 

• In the case of r as final lepton, we see sizable differences between all models' predic- 
tions. The nearest results to the SM correspond to the SUSY I and III. 
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BR{Af, — > Ae'*"e ) Region I Region II Region III 



SM 
SUSY I 
SUSY II 
SUSY III 
SUSY SO(IO) 



2.86 X 10-*5 1.12 X 10-6 0.81 x 10"^ 

4.57 X 10-6 1.56 X 10-6 0.99 x 10-^ 

5.66 x 10-6 2.69 x 10-6 2.61 x 10-6 

2.93 X 10-6 0.98 X 10^6 0.66 x 10-6 

2.65 X 10-6 1.20 X 10-6 1.09 x 10-6 



Table 6: The central values of branching ratio for — > Ae'^e decay channel at different 
regions in SM and different SUSY models. 



SM 



2.25 X 10-6 1.12 X 10-6 0.81 x 10-6 



SUSY II 4.65 X 10-6 2.63 x 10-6 2.55 x 10" 



SUSY SO(IO) 2.33 x 10-6 1.20 x 10-6 1.09 x 10-6 

Table 7: The central values of branching ratio for A^ A/i^/^- decay channel at different 
regions in SM and different SUSY models. 

3.3 The branching ratio 

Integrating the differential branching ratio over s in the considered regions and taking into 
account the central values of the form factors, we find the branching ratios for various 
models as presented in Tables 6, 7 and 8 for different lepton channels. A quick glance at 
these Tables leads to the following results: 

• as it is expected the values of the branching ratio decrease when going from the e to 

T. 

• The order of branching ratios indicates that these channels are accessible at the LHC. 
Note that as we have already mentioned this decay channel has been observed by 
CDF Collaboration at Fermilab in channel [7] . 

• All SUSY models have predictions considerably different than those of the SM in all 
regions and at all lepton channels. 

• The maximum deviation from the SM results belongs to the SUSY II model. When 
considering the numerical values, the maximum deviation of the SUSY II result from 
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BR{Ab At+t~) Region I Region II 



SM 


0.87 


X 


10- 


-7 


3.84 


X 


10- 


-7 


SUSY I 


1.35 


X 


10- 


-7 


5.55 


X 


10- 


-7 


SUSY II 


2.01 


X 


10- 


-7 


2.44 


X 


10- 


-6 


SUSY III 


0.72 


X 


10- 


-7 


2.13 


X 


10- 


-7 


SUSY SO(IO) 


1.12 


X 


10- 


-7 


1.61 


X 


10- 


-6 



Table 8: The central values of branching ratio for — > At'^t decay channel at different 
regions in SM and different SUSY models. 



the SM prediction corresponds to the region II for r channel. In this case, the result 
of the SUSY II is approximately 6 times greater than that of the SM. 



3.4 The FBA 

The lepton forward-backward asymmetry {Afb) is defined as: 

where Nf is the number of moving particles to forward direction and Nb is the number of 
moving particles to backward direction. In technique language, the lepton FBA is written 
in terms of the differential decay rate as: 

'z,s)dz— / {z,s)dz 



dsdz J_i dsdz 



^Ms) = '-^TT-^ '-it-^ ■ (3.22) 







dsdz 



r d^T 



Using this definition, we plot the dependence of the lepton FBA on s for e, ^ and r 
channels in the SM and different SUSY models in Figures 3 and 4. From these figures which 
are also plotted considering the central values of the form factors, it is clear that, 

• in the case of the e and ix channels, the SUSY I and II behave different than the 
other models. In these channels for the small values of the s, the maximum deviation 
belongs to the SUSY I, however, for higher values of the s the maximum deviation 
corresponds to the SUSY II. 

• In r channel, the maximum deviation from the SM prediction belongs to the SUSY 
III. 
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Figure 3: The dependence of the FBA on s for Af, — i- Ae^e and Af, — )■ A/i+/i transitions 
in SM and different SUSY scenarios using the central values of form factors. 
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Figure 4: The dependence of the FBA on s for A;, — t- Ar+r transition in SM and different 
SUSY scenarios using the central values of form factors. 



• The zero points of the FBA in different SUSY models move slightly to the left compare 
to the SM predictions. In some regions, the SUSY I, II and III have different signs 
with the SM predictions. 

• The SUSY SO(IO) represents overall the closest results to the SM predictions. 

3.5 The physical quantities under consideration taking into ac- 
count the uncertainties of the form factors 

In this subsection, we would like to consider the above mentioned physical quantities taking 
into account the uncertainties of the form factors and discuss the effects of these errors on 
the results. For this aim, considering the errors of the form factors, we plot the dependence 
of the differential branching ratio and FBA on s at different lepton channels and different 
models in figures 5-8. From these figures we see that. 
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Figure 5: Comparison of the results of differential branching ratio with respect to s for 
the Afo —7- Ae^e^ transition obtained from different SUSY models with that of the SM 
considering the errors of form factors. 
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Figure 6: The same as Figure 5 but for /x. 

• as far as the differential branching ratio are concerned, at the e and fi channels, the 
band of SUSY SO(IO) approximately covers the band of the SM. In the case of the 
SUSY I, II and III, although their bands coincide with that of the SM, there are some 
regions that these SUSY models have different predictions. Among these different 
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Figure 7: The same as Figure 5 but for r. 
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Figure 8: Comparison of the results of the FBA with respect to s for all lepton channels 
obtained from different SUSY models with that of the SM considering the errors of form 
factors. 



SUSY models, the maximum discrepancy from the SM prediction belongs to the 
SUSY II. In r channel, the difference between different SUSY models predictions and 
that of the SM can not be completely killed by the errors of form factors for any SUSY 
models. In this channel, there are also common regions between the bands of SUSY 
models and that of the SM, except than the SUSY II which does not approximately 
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nK[Ai, — >■ Ae^e j 


Region I 


T> TT 

Region II 


TJ TTT 

Region III 


SM 


(2.86 ± 1.43) X 10-6 


(1.12 ±0.56) X 10-6 


(0.81 ± 0.40) X 10-6 


SUSY I 


(4.57 ±2.42) X 10-6 


(1.56 ±0.82) X 10-6 


(0.99 ±0.52) X 10-6^ 


SUSY II 


(5.66 ± 3.05) X 10-6 


(2.69 ± 1.45) X 10-6 


(2.61 ± 1.40) X 10-6 


SUSY III 


(2.93 ± 1.52) X 10-6 


(0.98 ±0.50) X 10-6 


(0.66 ±0.34) X 10-6^ 


SUSY SO(IO) 


(2.65 ± 1.35) X 10-6 


(1.20 ±0.61) X 10-6 


(1.09 ±0.55) X 10-6 


Table 9: The values of branching ratio for A 
in SM and different SUSY models considerin 


6 — > Ae'^e decay channel at different regions 
g the uncertainties of form factors. 










SM 


(2.25 ± 1.12) X 10-6 


(1.12 ±0.56) X 10-6 


(0.81 ± 0.40) X 10-6 




SUSY II 


(4.65 ±2.51) X 10-6 


(2.63 ± 1.42) X 10-6 


(2.55 ± 1.37) X 10-6 




SUSY SO(IO) 


(2.33 ± 1.18) X 10-6 


(1.20 ±0.61) X 10-6 


(1.09 ±0.55) X 10-6 



Table 10: The values of branching ratio for A;, — )■ A/x"*"// decay channel at different regions 
in SM and different SUSY models considering the uncertainties of form factors. 



coincide anywhere with the SM result. 

• In the case of the FBA, the errors of form factors do not approximately affect the 
central values, except than the higher values of s, which we see narrow bands for 
different SUSY models as well as the SM. At e and fi channels, the bands of the 
SUSY I, SUSY III, SUSY SO(IO) and SM coincide with each other somewhere at 
higher values of s, but the SUSY II has different prediction. At r channel, all models 
have different predictions. 

• As it is expected, the forward-backward asymmetry and in particular its zero-crossing 
points are more robust than the differential branching ratio such that they are not 
approximately affected by the uncertainties of the form factors. This is the case also 
in B ^ K^*H+l- channel. 

Now, we discuss the effects of the uncertainties of the form factors for the branching 
ratios at different regions and for different models. Taking into account the errors of 

form factors, we present the vahics of branching ratios at different channels in Tables 
9, 10 and 11. From these Tables we deduce the same results as we have seen from the 
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nK[Ab — >■ Ar^T j 


T) T 

Region I 


Region 11 


SM 


(0.87 ±0.43) X 10-^ 


(3.84 ± 1.92) X 10-^ 


^ SUSY I 


(1.35 ±0.71) X 10-^ 


(5.55 ±2.94) X 10-^ 


CTTC!"V TT 


[2.[JL ± i.Uo j X lU 


± l.olj X lU 


SUSY III 


(0.72 ±0.37) X 10-^ 


(2.13 ± 1.10) X 10"^^ 


SUSY SO(IO) 


(1.12 ±0.57) X 10-^ 


(1.61 ±0.82) X 10-6 



Table 11: The values of branching ratio for — > Kt'^t decay channel at different regions 
in SM and different SUSY models considering the uncertainties of form factors. 

figures of the differential branching ratio. Although the central values for different 
SUSY models and the SM differ considerably from each other, considering the errors 
of the presented results in these Tables, we observe that approximately in all cases 
the results of different models coincide with each other, except for the SUSY II and 
SUSY SO(IO) at T channel and region II which have considerable discrepancy with 
the other model predictions. 

4 Conclusion 

In the present work, we have calculated the amplitude and differential decay rate for the 
semileptonic Af, — > A£+£- transition in different supersymmetric models. We have taken 
into account all twelve form factors entered the low energy matrix elements and recently 
calculated via light cone QCD sum rules in full theory to analyze the differential branch- 
ing ratio, total branching fraction and the lepton forward-backward asymmetry. We have 
considered different SUSY scenarios in the calculations and compared the obtained results 
with the SM predictions. As far as the central values of the form factors are considered, 
in general, we observed considerable deviations from the SM predictions. In the case of 
the (differential) branching ratio, the maximum deviations from the SM predictions belong 
to the SUSY II scenario. As far as the FBA is concerned, at e and ji channels and lower 
values of the s, the maximum deviation belongs to the SUSY I, however, for the higher 
values of the s and the same lepton channels, the maximum discrepancy corresponds again 
to the SUSY II model. Taking into account the uncertainties of the form factors, we have 
observed that the branching ratio is more affected by these errors. The bands of the SUSY 
SO(IO) approximately cover the SM bands at the e and ji channels. For other SUSY models 
and all lepton channels, although we have seen some intersection regions between different 
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SUSY bands and the SM predictions, there are considerable discrepancies between the SM 
and SUSY models predictions. Especially, at r channel, there is a big discrepancy between 
the SUSY II and the SM bands. When we consider the FBA, the uncertainties of the form 
factors do not affect this quantity and its zero-crossing points. We see overall a consid- 
erable discrepancies between the narrow bands of the different considered SUSY models 
and that of the SM. Such discrepancies can be considered as a signal for existence of the 
supersymmetric particles. 

The orders of the branching ratio at all lepton channels and all the considered regions 
of depict that these decay channels can be checked at LHC in near future. Note that 
as we have also previously stressed, this channel for \x case has been observed recently by 
CDF Collaboration at Fermilab. We are waiting for the LHCb Collaboration results on 
these channels, which they have in their physics program [8]. 

Comparison of the experimental results on the branching ratio as well as the FBA with 
the predictions of the present work, especially determination of the sign and zero-crossing 
points of the FBA, which have not been affected by the errors of the form factors, can help 
us get valuable information about the existence of the SUSY particles. 

As we have already noticed, in numerical analysis, we have used the values of the Wilson 
coefficients presented in Table 5 for different SUSY scenarios. These values are obtained 
when the masses of the neutral Higgs bosons are taken in the interval (91 — 200) GeV 
(see for instance [4, 15-17]). Considering the recent developments by the CMS and ATLAS 
Collaborations at CERN on the mass of the Higgs-like boson (~ 125 GeV), the used values 
of the Wilson coefficients are still viable. However, after clarifying whether the obtained 
boson at LHC is the standard or non-standard Higgs, and using its exact mass, one may 
recalculate the Wilson coefficients in different SUSY scenarios. Obviously, it will be possible 
to improve the obtained results in the present work using the new values of the Wilson 
coefficients. 

References 

[1] J. Incandela, CMS talk at "Latest update in the search for the Higgs boson at CERN" , 
July 4, 2012; F. Gianotti, ATLAS talk at "Latest update in the search for the Higgs 
boson at CERN", July 4, 2012. 

[2] M. Hirsch, F. R. Joaquim, A. Vicente, "Constrained SUSY seesaws with a 125 GeV 
Higgs" [ arXiv: 1207.6635 [hep-ph]]. 



20 



[3] Q.-Sh. Yan, Ch.-Sh. Huang, L. Wei, Sh.-H. Zhu, "Exclusive Semileptonic Rare Decays 
B {K,K*)e+e- in Supersymmetric Theories", Phys. Rev. D 62, 094023 (2000) 
[arXiv:hep-ph/0004262]. 

[4] W.-J. Li, Y.-B. Dai, Ch.-Sh. Huang, "Exclusive Semileptonic Rare Decays B — )■ 
K^*H+l- in a SUSY SO(IO) GUT", Eur. Phys. J. C 40, 565 (2005) [arXiv:hep- 
ph/0410317]. 

[5] M. J. Aslam, Y.-M. Wang, C.-D. Lii, "Exclusive semileptonic decays of A;, — >■ AZ+Z~ in 
supersymmetric theories", Phys. Rev. D 78, 114032 (2008) [arXiv:0808.2113 [hep-ph]]. 

[6] A. Ahmed, I. Ahmed, M. A. Paracha, M. Junaid, A. Rehman, M. J. Aslam, "Compar- 
ative Study of Be — > D*£'^i~ Decays in Standard Model and Supersymmetric Models" 
[ arXiv:1108.1058 [hep-ph]]. 

[7] T. Aaltonen et al. [CDF Collaboration], "Observation of the Baryonic Flavor-Changing 
Neutral Current Decay Ag A/^V"", Phys. Rev. Lett. 107, 201802 (2011) 
[arXiv:1107.3753 [hep-ex]]. 

[8] Our personal communications with Yasmine Sara Amhis from LHCb Collaboration. 

[9] J. Beringer et al. (Particle Data Group), Phys. Rev. D 86, 010001 (2012). 

[10] T. M. Aliev, K. Azizi, M. Savci, "Analysis of the A,, A£+£- decay in QCD", Phys. 
Rev. D 81, 056006 (2010) [arXiv:1001.0227 [hep-ph]]. 

[11] G. Buchalla, A. J. Buras, M. E. Lautenbacher, "Weak decays beyond leading loga- 
rithms". Rev. Mod. Phys. 68, 1125 (1996) [arXiv:hep-ph/9512380]. 

[12] M. Beneke, G. Buchalla, M. Neubert, C. T. Sachrajda, "Penguins with charm and 
quarkhadron duality", Eur. Phys. J. C 61, 439 (2009) [arXiv:0902.4446 [hep-ph]]. 

[13] A. Khodjamirian, Th. Mannel, A.A. Pivovarov, Y.-M. Wang, "Charm-loop effect in 
B K^*k+£- and B K*i\ JHEP 1009, 089 (2010) [arXiv:1006.4945 [hep-ph]]. 

[14] M. J. Aslam, C.-D. Lu, Y.-M. Wang, "B K^{U30)l+l- decays in supersymmetric 
theories", Phys. Rev. D 79, 074007 (2009) [arXiv:0902.0432 [hep-ph]]. 

[15] J.-F. Cheng, Ch.-Sh. Huang, X.-H. Wu, "Neutral Higgs boson contributions to CP 
asymmetry of £ ^ ^Ks in MSSM", Phys. Lett. B 585, 287 (2004) [arXiv:hep- 
ph/0306086]. 



21 



[16] Ch.-Sh. Huang, P. Ko, X.-H. Wu, Y.-D. Yang, "MSSM Anatomy of the Polarization 
Puzzle in S ^ (j)K* Decays", Phys. Rev. D 73, 034026 (2006) [arXiv:liep-pli/0511129]. 

[17] J.-F. Cheng, Ch.-Sh. Huang, X.-H. Wu, "CP asymmetries in B ^ (j)Ks and B r] Kg 
in MSSM", Nucl. Phys. B 701, 54 (2004) [arXiv:hep-ph/0404055]. 



22 



